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We consider the self-localization of neutral impurity atoms in a Bose-Einstein condensate in a 
ID model. Within the strong coupling approach, we show that the self-localized state exhibits 
parametric soliton behavior. The corresponding stationary states are analogous to the solitons of 
non-linear optics and to the solitonic solutions of the Schrodinger-Newton equation (which appears 
in models that consider the connection between quantum mechanics and gravitation). In addition, 
we present a Bogoliubov-de Gennes formalism to describe the quantum fluctuations around the 
product state of the strong coupling description. Our fluctuation calculations yield the excitation 
spectrum and reveal considerable corrections to the strong coupling description. The knowledge of 
the spectrum allows a spectroscopic detection of the impurity self-localization phenomenon. 
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I. INTRODUCTION 

Recent work pointed out that neutral impurity atoms 
immersed in a dilute gas Bose-Einstein condensate (BEC) 
can spontaneously self- localize 0, . Repeated measure- 
ments of the impurity position would observe this par- 
ticle's wavefunction to have a finite extent, even if its 
average position cannot be predicted a priori (at least 
if the BEC were homogeneous in the absence of the im- 
purity). In three dimensions (3D), the BEC-impurity 
physics is akin to that of an electron moving in a po- 
lar crystal (a polaron) Q. In 3D, the self-localization 
of a neutral atom BEC-impurity occurs when the mag- 
nitude of the impurity-boson scattering length exceeds 
a minimal value (regardless of the sign of that interac- 
tion). The self- localized BEC impurity state resembles 
that of a small polaron. Traditional (i.e. electronic) small 
polarons have been described successfully in the strong 
coupling limit by the Landau-Pekar treatment j^J which 
assumes a wavefunction that is the product of a single 
electron wavefunction and a bosonic phonon state. The 
self-localized BEC-impurity has been described similarly 
[1,2] by the product of a single particle impurity wave- 
fimction and a BEC-state. Linearizing the BEC-response 
to the impurity density then gives an effective impurity 
equation of a single particle experiencing an attractive 
Yukawa self-interaction potential [1]. While the ensuing 
analysis is elegant, the validity of the starting point (i.e. 
the product state) is much more questionable than in the 
description of traditional polarons. In the latter case, the 
large mass difference of crystal ions and electrons implies 
a clear separation of time scales which justifies the prod- 
uct state - one can always think of the slow (boson) field 
adjusting itself to the time averaged field of the fast (elec- 
tron) particle. In the cold atom BEC-impurity, the boson 
and impurity atoms tend to have similar masses and the 
time scales do not separate. Hence, we can expect the 
deviations from the product state, the fiuctuations that 



describe the entanglement of the impurity and boson de- 
grees of freedom, to become much more significant. 

In this paper, we describe the one-dimensional ana- 
logue of the BEC-impurity, realizable in quasi-one- 
dimensional BEC's. For this system and for a specific 
choice of the parameters, we find an exact solution to the 
strong coupling equations and we solve for the quantum 
fluctuations as a function of the impurity-boson mass ra- 
tios. The explicit solution to the strong-coupling equa- 
tions are parametric solitons, which establishes that the 
ID analogues of self-localized BEC- impurities are optical 
solitons. Specifically, the impurity soliton solutions are 
resemble the solitons that appear when a quadratic non- 
linearity is responsible for second-harmonic generation 
1^1^. The Schrodinger-Newton equation that models the 
gravitational interaction while preventing the quantum- 
mechanical spreadingof the center-of-mass position of 
macroscopic objects S ; as well as the mean-field 
equations of coupled atomic- molecular BEC's 0, ^ 
possess similar parametric soliton solutions. The con- 
nection with these very different fields of physics further 
extends the scope of the BEC-impurity physics. 

The experimental realization of the BEC-impurity sys- 
tems requires the creation of distinguishable atoms in 
BEC's. This feat has been realized by converting a 
fraction of the BEC-atoms with a two-photon Raman 
transition to a different spin state in order to observe 
the superfluid suppression of slow impurity scattering by 
the BEC or by trapping distinct species of atoms 
[TsL IT^ or isotopes jl^]. In this paper, we describe a 
BEC trapped in a one dimensional box. Our predictions 
apply to atomic traps with strong confinement in the two 
transverse directions (quasi- ID). A ID box potential was 
recently achieved experimentally but our calcula- 
tions should also describe the physics of quasi ID-BEC's 
with trapping potentials that vary slowly in the longitu- 
dinal direction. 

The paper is organized as follows: In section II, we 
introduce the model and work within the product state 
description. Section III presents a formalism that goes 
beyond the product state ansatz by means of a Bogoli- 
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ubov description of the quantum fluctuations. In section 
IV, we discuss our numerical results, and we conclude in 
section V. 



II. STRONG COUPLING APPROACH 

We consider M impurity bosonic atoms immersed in a 
homogeneous BEC in a ID model. The Hamiltonian of 

the system reads 
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where (t){x) and V'(a^) stand for the condensate boson and 
impurity atom field operators, respectively; ttib and nii 
are the boson and impurity masses, while /iB and Ei rep- 
resent the chemical potentials of the BEC and impurity 
systems. We assume contact boson-boson and impurity- 
boson interactions characterized by Abb and Aib , respec- 
tively. We neglect the mutual interactions of impurity 
atoms in the assumption that their number and local 
density remains sufficiently small. 

In the strong coupling treatment we describe the 

system in terms of a product state. 
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where 77(x) is the wave- function occupied by the M impu- 
rity atoms and ({x) is the single particle state occupied 
by the N bosons. The expectation value of the Hamil- 
tonian ^ for the product state Q gives an energy that 
reaches its minimal value when the following equations 
are satisfied: 
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and we have assumed that both ip{x) and V'o(2;) are real- 
valued - an assumption that is permitted since we wish 
to describe a ground state. We also assume that M, the 
number of impurity atoms is small, so that we may ex- 
pect only a slight modification of the boson wave- function 
ip{x) with respect to the homogeneous BEC solution 
where p is the BEC density. Therefore we substitute 



'fix) = Vp + <^o(a;), 
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into Q and keep linear terms in ipoix) only, 



2mB 



2AbbP 



-—^ + 2XiB^(t)o{x 
2mi 



Mx) + Xm^P'>Pl{x)^Q, (6) 
Mx)^Ei4,o{x), (7) 



where /iB ~ AbbP and Ei = Ei — XmP !]• This approxi- 
mation is valid provided. 
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The Eq. (O can be solved in terms of the Green function 
of the ID Hclmholtz equation. 
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represents the BEC-coherence length. The substitution 
of © into Eq. lO gives an equation that describes im- 
purity atoms self-interacting through an attractive expo- 
nential potential. Note the difference with the 3D situ- 
ation in which the self-interaction takes place through a 
Yukawa (screened Coulomb) potential 1]. 

Equations similar to Eqs. (|6I7|) occur in non-linear op- 
tics 0, , in mean-field descriptions of coupled atomic- 
molecular BEC's 0,^3 and in the Schrodinger-Newton 
model 0. These equations are known to possess para- 
metric soliton solutions. In the present case, an analyti- 
cal solution exists for a particular value of the condensate 
density. For 
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These coupled equations possess a solitonic solution, 
3 ' 
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In Eqs. (|12I13|I we have scaled energy and length by Eq 
and xq where 
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The soliton H13|l describes a bound state of the impurity- 
BEC system with an impurity extend of order y^6/M, 
which is precisely the BEC healing length since 
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The self-localization takes place for either sign of the 
impurity-BEC interaction (Aib > and Aib < 0). In 
addition to and caused by the impurity localization, the 
condensate is deformed. The BEC exhibits a dip in the 
density if the impurity and BEC-particle mutually repel 
(Aib > 0) and a hump if the they attract (Aib < 0). 

We have shown there exists a class of analytical so- 
lutions for a particular choice of parameters of the ID 
BEC-impurity system. In general, one can solve the cou- 
pled equations numerically. Using the Gaussian ansatz 
for "00(2;), one can also show that there is impurity self- 
localization in the ID-model even as Aib -> [13. This 
is markedly different from the 3D-situation for which the 
impurity-BEC interaction must be sufficiently strong be- 
fore self-localization sets in P|. 



then caused entirely by the interactions with the impuri- 
ties. We expect that this approximation will not greatly 
affect the impurity physics. It does, however, modify the 
description of long wavelength BEC-excitations. In the 
absence of impurities the approximation H17() corresponds 
to the Hartree-Fock approximation. This approach pre- 
dicts a gap in the BEC-excitation spectrum, contrary to 
the Goldstone (or Hugenholtz-Pines) theorem, and dif- 
ferent from the the Bogoliubov spectrum ■ The final 
effective Hamiltonian in the units (|14ll reads 
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We treat the resulting Hamiltonian (|18() with a Bogoli- 
ubov approximation [ig: first we solve the mean-field 
equations, then we construct the quantum fluctuations 
around the mean-field solutions. To this end, we expand 
the field operators as 



III. QUANTUM FLUCTUATIONS 

Hamiltonian ^ describes a small number of bosonic 
impurity particles immersed in a BEC in a ID-model 
with inter-particle interactions of the contact type. As- 
suming that the presence of the impurity atoms does not 
perturb the condensate significantly, we may decompose 
the bosonic field operator into 



if{x) « ^ + 4>{x), 



(16) 



where denotes the stationary mean-field solution for 
a condensate wavefunction in the absence of impurities, 
and where (t>{x) describes the small perturbations of the 
condensate caused by the impurity (or impurities). We 
substitute IjlGI) into the Hamiltonian and keep terms 
that are quadratic in the operator only. Next, we re- 
place the condensate boson-boson interaction term in the 
Hamiltonian by 
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This approximation implies that we neglect the depletion 
of the BEC induced by the interactions between conden- 
sate particles. The BEC-depletion that we compute is 



ip{x) 
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In zeroth order in the S(j) and dip operators one obtains 
the mean field equations. 
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identical to Eqs. Ht)l7|) . written in the units (|14|l . 

While the mean-field description of coupled atomic- 
molecular BEC's give solitonic solutions that are identi- 
cal to Eq. (|13|l |1(1 , the Hamiltonian of the atomic- 
molecular BEC's is different from Eq. H18|l . The coupled 
atomic-molecular BEC Hamiltonian contains terms that 
convert molecules into atoms and vice versa, rather than 

the interaction term 7 (p'' (x) + (t){x) il;'^{x)ip{x). The 
latter term can be found in models where bosonic parti- 
cles ip feel a long-range force caused by the exchange of 
a mesonlike particle (f) [TH . 

Given that "00 (2;) and (poi^) satisfy the mean-field 
equations H21|l . the first order terms in the d(j) and Sip 
operators cancel. The second order term gives an effec- 
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tive Hamiltonian that takes the form 
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+ j{Mx)+roi^)}-E. (24) 



The problem of diagonahzing the Hamiltonian (|22|l re- 
duces to the problem of diagonahzing the non-Hermitian 
operator £ (i.e. to solving equations of the Bogoliubov- 
de Gennes type). The ^-operator possesses two sym- 
metries (similar to the symmetries of the original 
Bogoliubov-de Gennes equations 0), 
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and 
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(26) 



(27) 



are the first and third Pauli matrices, respectively. 

Suppose that all eigenvalues of the /^-operator are real. 
The symmetries ()25|l imply that if 
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V \vt) J 
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is a right eigenvector of the >C-operator with eigenvalue 
£k, then l^*^) = wsl^jl?) is a left eigenvector of the same 
eigenvalue e^, and uil^f*) is a right eigenvector with 
eigenvalue — e^. Except for the eigenstates corresponding 
to zero eigenvalue, the eigenstates of the >C-operator can 
be divided into two families and "— ", 
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We apply the Bogoliubov transformation 
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with quasi-particle operators that can be written as 



bk = (ut\<j> 



\^) + {ut\^^)-{vt\^) + {ut\^^). 

-K^I<^^> + (4l^>-Kl^^), (31) 



and that fulfill the bosonic commutation relation 
= (5fc,fc'- This transformation gives an effective 
Hamiltonian that is diagonal 



(32) 
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To obtain the energy eigenvalues and eigenstates of the 
effective Hamiltonian (|22|) . one has to solve the mean- 
field equations 1)21(1 . then diagonalize the operator H23|l. 
For a specific value of the average BEC-density, given 
by Eq. Hll|l . we obtain the solitonic solution ((13|l of the 
mean-field equations. In that case, the eigenstates of the 
£-operator corresponding to zero-eigenvalue take on the 
form 
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- cosh 







-2 



dx cosh \ — 
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(33) 



where the first eigenvector corresponds to the break- 
ing of the U(l) symmetry in the BEC-Bogoliubov the- 
ory [Tflj l , while the second eigenvector corresponds to the 
breaking of the translational symmetry, indicating that 
the translation of the soliton costs no energy. There is 
another zero-momentum eigenstate which has non-zero 
eigenvalue. This eigenvalue and its eigenvector take on 
simple analytical forms. 



Aim 




,(34) 



where L is the size of the ID-box. Other eigenstates can 
be found numerically. In the next section, we present 
numerical results for realistic parameter values. 
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FIG. 1: Excitation spectrum of the quasi ID Na BEG with 
eight *^Rb atom impurities, in the energy units of Eq. 1141 . 
The upper branch corresponds to BEG-particle like excita- 
tions, while the impurities remain localized. This branch co- 
incides with the Hartree-Fock spectrum of the BEG without 
the impurity atoms (diamonds) given by Eq. The lower 

branch corresponds to excitations that delocalize the impu- 
rity atoms (except for the lowest energy level, as can be seen 
in Fig|5J. Gircles denote even solutions and crosses show so- 
lutions of odd symmetry. 



FIG. 2: The eigenstates of the /^-operator, Eq. 12311 . corre- 
sponding to the lower energy branch of the Hamiltonian 1321 , 
see Fig. Q Left column is related to condensate modes, i.e. 
uf. (solid lines) and v'^ (dashed lines — hardly visible behind 
the solid lines), while the right column is related to impuri- 
ties modes, i.e. uf (solid lines) and (dashed lines). There 
are four lowest energy eigenstates presented in the figure, i.e. 
(a)-(b) is the lowest one, (c)-(d), (e)-(f) and (g)-(h) are shown 
in order of increasing energy. The states describe the excita- 
tion of eight self-localized *^Rb atoms embedded in a quasi 
ID ^^'Na BEG. 



IV. NUMERICAL RESULTS 

We consider eight impurity *^Rb atoms, i.e. M — 8, 
immersed in a BEC of ^^Na atoms in an elongated trap 
that is a box in the axial direction while in the trans- 
verse directions there is a harmonic trap of frequency 
Lu±. The transverse ^^-confinement is so strong that 
only the ground states of the transverse degrees of free- 
dom are relevant. In this quasi ID-regime, the system 
is effectively one dimensional and confined by a 1 D-box 
potential. The coupling constants of the Hamiltonian 
read 

i 9 fi2 / 1 1 \ 1 

AiB ~ 27rft am 1 , n , — 2T 

\mi TUB J 7r(CTg + erf) 

= 2hujj_aiB, 

Abb = 7: — 2" = 2?iwxaBB, (35) 

TTlB 27rCTg 

where the s-wave scattering lengths obb = 3.4 nm 
and aiB = 16.7 nm The ctb = \JliiTrvBi^ and 

a\ = fi / miuj ± lengths represent the ground state ex- 
tents of a single BEC-boson particle and of a single im- 
purity atom confined by the two-dimensional harmonic 
trap in the transverse direction with frequency We 
take the transverse trapping frequency to be equal to 
uj±_ = 2tt X 500 Hz. Assuming that the size of the box 
in the longitudinal direction is equal to lOO/xm (in units 
(|14|l . this length corresponds to a box size L « 80), we 
confine 11000 atoms in the BEC. These parameters cor- 
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FIG. 3: The same as in Fig. |5| but for the upper energy 
branch of the Hamiltonian 1321 . see Fig.Q The modes shown 
in panels (a)-(b) correspond to the solution given by 1341 . 



respond to the specific value of the condensate density 
(|ll|l and are well within the experimental range . 

Figure Q shows the excitation spectra that were cal- 
culated numerically with periodic boundary conditions. 
Note that the dispersion shows two branches. Each exci- 
tation of the lower branch delocalizes an impurity atom, 
except for the very lowest level, which corresponds to 
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FIG. 4: Panel (a): solid line denotes the density of the 
impurity atoms, iPq{x), calculated within the product state 
ansatz, dashed line is the correction to this density, i.e. 

{0\5'ilj^x)Sip{x)\0) = J2k bti^)]^^ obtained within the for- 
malism that goes beyond the product state approximation. 
Panel (b): BEG density given by Eq. (I39II . The results are re- 
lated to eight *^Rb atoms immersed in a BEG of 2^Na atoms. 
The length is given in the units 1141 . 



an excitation that leaves the impurity atom localized, 
as shown in Fig. [5] The upper branch corresponds to 
particle-like BEC excitations in which the impurity re- 
mains localized, as shown in Fig. O The highest curves 
of Fig.Oshow the mode of Eq. (|34|) . In Fig.^we have also 
plotted the BEC-spectrum calculated within the Hartree- 
Fock approximation p^ . 
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where n — 0, 1, . . .. Actually, to obtain H36|l we have 
used the Hamiltonian (|18|l in which we neglected terms 
involving the ip{x) operator. The Hartree-Fock spec- 
trum H36|l provides the relevant comparison for the model 
since the terms neglected in Ea. (|17|) correspond to the 
Hartree-Fock description of the single BEC. Figure 
shows that the upper energy branch is nearly identical to 
the Hartree-Fock BEC-spectrum in the absence of impu- 
rities. The lower energy branch describes impurity exci- 
tations with negligible excitation of the BEC, see Fig.lSJ. 
We interpret the energy gap in the lower energy branch 
(as the wavenumber tends to zero) as the minimal energy 
needed to break the self-localization bond of the many- 
body BEC-impurity system. This gap is then a char- 
acteristic feature of impurity self-localization. Its detec- 
tion by means of Bragg spectroscopy j23|, for instance, 
can serve as a smoking gun signal for the observation of 
BEC-impurity self-localization. 

The ground state of the system is the Bogoliubov vac- 
uum state |0), i.e. the state annihilated by all bk op- 
erators We can estimate the quantum fluctuation 
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FIG. 5: The integral (I38II [panel(a)] and condensate deple- 
tion [panel (b)] as a function of the impurity mass mi. Note 
that the maximal values of the and condensate depletion 
correspond to mi/rriB ~ 0.5. 



corrections to the product state by calculating the den- 
sity of the impurity atoms. 



{0\ip^x)'ip{x)\0) = t/jlix) + {0\5'ip''{x)SiP{x)\0) 
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The quantum 



In Fig. ^ we plot ijQ{x) and J2k 
fluctuation contribution to the impurity density is 

2 
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= 0.31, 



(38) 



as compared to the value 8 for the integral of ip'^{x), 
which shows that the fluctuation corrections to the prod- 
uct state are not insignificant. In the case of a single im- 
purity, i.e. M — 1^ keeping the other parameters fixed, 
C^p — 0.30 and the fluctuations corrections turn out to 
be more dramatic. The condensate density (see Fig. 0} 
reads, 



(0|^t(^)^(^)|0) = 



[^/x^ + (t>n{x)] 
{Q\6^\x)64>{x)\{)) 
xop + 2y^x^(j)o{x), 



(39) 



where we have introduced approximations due to the 
facts that (1)q{x) / ^x^p < 0.083 [that justifies also the 

linearization © or (O] and (O|50t(o)50(o)|O)/(/)§(O) w 
0.03. 

Decreasing the mass of the impurity atom, while keep- 
ing all other parameters fixed, the integral and the 
BEC depletion increase and reach their maximal values 
for mi/wB « 0.5, see Fig. |S1 This shows that when the 
impurity mass is half that of the BEC-bosons, the prod- 
uct state approximation needs the largest corrections. 
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Qualitatively, this behavior agrees with the time scale 
separation argument, which predicts that the large and 
small mass ratio regimes are well described by a product 
state 

V. CONCLUSIONS 

We have considered the self-localization of neutral, 
bosonic impurity atoms embedded in a dilute gas Bose- 
Einstein condensate in a ID-model. We have analyzed 
the system within a Bogoliubov formalism that describes 
the quantum fluctuations around the strong coupling 
product state approximation previously made in the cold 
atom literature QQ]- Our description gives an excitation 
spectrum that consists of two branches. The lower energy 
branch corresponds to excitations that delocalize the im- 
purity atoms (with the exception of the lowest energy ex- 
citation). The higher energy branch corresponds to BEC- 
particle like excitations (and coincides, in our approxima- 
tion, with the Hartree-Fock spectrum of the BEG in the 
absence of the impurity atoms) while the impurities re- 
main self-localized. The energy gap in the lower energy 
branch (i.e. the branch that corresponds to the excitation 
of impurity atoms with negligible excitation of the BEG) 
suggests a spectroscopic means (Bragg spectroscopy pl|) 
for detecting BEG-impurity self-localization. 

The parameters of our calculation are well within the 



experimental range. An extension of this approach to the 
three-dimensional BEC-impurity system appear straight- 
forward. One can extend the formalism to account for the 
presence of a larger number of impurity atoms by includ- 
ing a nonlinear term in the Hamiltonian |^ to account 
for interactions among the impurity atoms. Also the de- 
pletion of the BEC caused by the interactions among con- 
densate particles (neglected in the present calculations) 
can be included. 

Furthermore, our calculations reveal the close anal- 
ogy of ID BEC-impurity self-localization with paramet- 
ric solitons known in nonlinear optics |^ in the 
Schrodinger-Newton model 0, 13, B and in the coupled 
atomic-molecular BEG svstem |lOl [TTI |. The problem 
of impurity atoms immersed in an atomic BEG offers 
intriguing opportunities for experimentally realizing the 
phenomena predicted in these fields. 
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